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Abstract 

Extensive numerical analysis of the eigenspectra of the SUq{N) in- 
variant Perk-Schultz Hamiltonian shows some simple regularities for 
a significant part of the eigenspectrum. Inspired by those results we 
have found two set of solutions of the associated nested Bethe-ansatz 
equations. The first set is obtained at a special value of the anisotropy 
(q = exp{i2-ir{N — 1)/N)) and describes in particular the ground state 
and nearby excitations as a sum of free-fermion quasienergies. The sec- 
ond set of solutions provides the energies in the sectors whose number 
Ui of particles of distinct species {i = 0, . . . , N — 1) are less or equal 
to the unity except for one of the species. For this last set we obtain 
the eigenspectra of a free fermion model for arbitrary values of the 
anisotropy. 

1. Introduction 

The nested Bethe ansatz is probably the most sophisticated algebraic 
construction of eigenvectors for integrable lattice models, where the underly- 
ing quantum group is of rank larger than one. In the same way as the usual 
Bethe ansatz, which was widely used during more then 60 years (see, e. g.. 
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|T]] for reviews) the amplitudes of the wave functions are expressed in terms 
of a sum of plane waves whose wave numbers are given in terms of non-linear 
and highly non trivial coupled equations known as the nested Bethe-ansatz 
equations (NBAE). 

Except in the thermodynamic limit, in very few cases the associated 
NBAE can be solved analytically. The exceptions happen for a small number 
of "spin waves" or for very small lattices. In addition to these, in the case 
where the equations are of non-nested type, as in the XXZ chain, analytical 
solutions for chains of arbitrary length can also be derived for the free fermion 
point (anisotropy A = 0) and for the special anisotropy A = — ^ ^ |[. 

In a recent paper we presented a new set of analytical solutions of 
the NBAE for finite quantum chains. These announced solutions correspond 
to the SUg{3) Perk-Schultz model 0, or the anisotropic SU{3) Sutherland 
model [0], with periodic or open boundary conditions, at a special value 
of the anisotropy. These solutions are given in terms of two types of wave 
numbers, however their number (po,Pi) do not satisfy the usual bounds {po < 
Pi < L), and it is not clear if although being solutions of the NBAE they 
do not correspond to zero-norm states. Moreover, even in the case where 
they correspond to physical eigenstates of the quantum Hamiltonian it is not 
clear, in general, in what eigensector should we expect these solutions. In 
the case of periodic boundary conditions, by using the functional equations 
of the model we were able to answer partially the above points in [Q. 

In the present paper we fill the above deficiency by presenting a set of 
new analytical solutions of the NBAE having the number of wave numbers 
{po,Pi) satisfying the usual bounds po < pi < L, and consequently explaining 
the free-fermion part of the eigenspectrum of the SU (3) Perk-Schultz model 
conjectured from extensive numerical calculations in ^j. Moreover our fur- 
ther "experimental" (i. e. numerical) investigations lead to new observations 
concerning the anisotropic SU{N) Perk-Schultz quantum chains. Based on 
our numerical evidences we again have formulated several conjectures. Sur- 
prisingly, the inspirations coming from these numerical observations enable 
us to find analytically a lot of solutions of the related NBAE, that explains 
part of these new conjectures. The list of analytical solutions we described 
in this paper are: 

• the sectors [riQ, rii, ^2) = {k,k + l,L — 2k — l) of the SU{3) model with 
q = exp(2z7r/3) with periodic boundary conditions. 
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• the sectors {k, k,L — 2k) and {k, k + 1, L — 2k — 1) of the SU{3) model 
with q = exp(2i7r/3) with free boundary conditions. 

• the sector (wq, rii, . . . , n^^i) = (1, 1, . . . , 1, L — + 1) of the SU{N) 
model with arbitrary values of q and free boundary conditions. 

• the sector {uq, rii, . . . , n7v-i) = (1, 1, . . . , 1, 2, L — A^) of the SU (N) 
model with q = exp((A^ — l)m/N) and free boundary conditions. 

• the sector (no, ni, n2, na) = (1,2,2,L — 5) of the SU{A) model with 
q = exp(?37r/4) and free boundary conditions Q 

The paper is organized as follows. In §2 we introduce the model, and give 
the NBAE for periodic and free boundary conditions. In §3 we consider the 
set of analytical solutions related to the SU{3) model with q = exp(2i7r/3), 
with periodic and free boundary conditions. In §4 we consider the NBAE for 
the sector (1,...,1,L — + with arbitrary values of the anisotropy. We 
consider in details the SU (4) case and formulate a mathematical statement 
that allows to generalize our consideration to the case of arbitrary rank. 
Using a special case of this statement in §5 we find the NBAE free-fermion 
solutions for the SU{N) model with anisotropy q = exp((iV — l)m/N) in 
the eigensector (1, . . . , 1, 2, L — A^). §6 contains a set of new "experimental" 
conjectures that merged from extensive numerical diagonalizations of the 
SU{N) quantum chains for arbitrary values of N. Finally in §7 we present 
our conclusions and a summary of our results. We also present in appendix 
A and B the solutions in the sector (1, 2,2, L — 5) of the S'f/(4) model with 
anisotropy q = exp(z37r/4) and the wave functions related with the solutions 
considered in §4, respectively. 

2. The SU{N) Perk-Schultz model 

The SU (N) Perk-Schultz model is the anisotropic version of the SU{N) 
Sutherland model with the Hamiltonian, in a L-site chain, given by 

L-l 

Hg = T.Hj,j+i+pHL,i (p = 0,l), where (1) 

N-l N-l 

Hid = - E E {EfEf + El^'Ef - qE^Ef - l/g^f ). 

a=0 b=a+l 

^We suspect that this class of solutions can be greatly generalized for the SU [N] model 
with q = exp((Af — \)iTT/N) (see discussion in the last section). 
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The N X N matrices E"''' have elements {E°-^)f.d = ^c^'d ^^"^ 1 ~ exp(i?7) 
plays the role of the anisotropy of the model. The cases of free and peri- 
odic boundary conditions are obtained by setting p = and p = 1 in (|I]), 
respectively. This Hamiltonian describe the dynamics of a system containing 

classes of particles (0, 1, — 1) with on-site hard-core exclusion. The 
number of particles belonging to each specie is conserved separately. Conse- 
quently the Hilbert space can be splitted into block disjoint sectors labeled 
by (no, ni, niv_i), where = 0, 1, L is the number of particles of specie 
i (i=0,l,...,N-l). The Hamiltonian (|l|) has a Sn symmetry due to its invari- 
ance under the permutation of distinct particles species, that implies that 
all the energies can be obtained from the sectors (ng, ni, n^v-i), where 
no < ... < riM-i and no + ni + ... + n7v-i = L. 

At g = 1 the model is SU{N) invariant and for q 1 the model has a 
U{1)^^^^ symmetry as a consequence of the above mentioned conservation. 
Moreover in the special case of free boundaries [p = 0), the quantum chain 
([l|) has a larger symmetry, being invariant under the additional quantum 
SU{N)q symmetry. This last invariance implies that all the eigenenergies 
belonging to the sector {nQ,n[, ...,n'j^_^) with nj, < n[ < ... < n'^_^ are 
degenerated with the energies belonging to the sectors (no, ni, njv-i) with 
no < Til < ... < niv-i, if Uq < tiq and Uq + n[ < no + rii and so on up to 
nQ + n[ + . . . + n^_2 < '^o + '^i + ■■■nN-2- 

The NBAE that give the eigenenergies of the SUq{N) Perk-Schultz model 
in the sector whose number of particles is (n^, i = 0, . . . , N — 1) are given by 
(seee. g. §,1) 

Pk Pk-l Pk+1 

n F(«fuf ) = n /(^fur^^) n /(^f\^r'^), (2) 

where k = 0,1, ...,N — 2 and i = l,2...,pk. The integer parameters pk are 
given by 

k 

Pk = ^ni, k = 0,l,...,N -2] p_i = 0, pn-i = L. (3) 

1=0 

The functions F{x,y) and f{x,y) are defined by: 

sin(3: -y-r]) sm{x - y -r]/2) 

F(x,y) = -, f(x,y) = j-r, (4) 

sm(x — y + T]) sm[x — y + ri/2) 
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for periodic boundary conditions, and by 



cos(2y) - cos(2x - 2r]) cos(2y) - cos(2x - r]) 

F{x,y) = — -— — — — , f{x,y) = — -— — ■ — -, 5 

cos(2?/) — cos(2x + 2ri) cos(2y) — cos(2x + i]) 

for the free boundary case. In using the NBAE (Q) we deal with variables 
of different order {wj^''} (fc = 0, . . . , iV — 2). The number of variables uf'^ of 
order k is equal to pk- The whole set of NBAE consists of subsets of order k 
which contain precisely pk equations {k = 0,1, N — 2). 

The eigenenergies of the Hamiltonian (|I|) in the sector (no, ni, n^v-i) 
are given in terms of the roots {"Uj^ ^^}: 

h\ q sm{uj+ri/2) sm{uj - ri/2)J ' 



where to simplify the notation we write uj = u^f^ (j = 1, . . . ,pAr_2). 

All the solutions of the NBAE (D which are going to be described in this 
paper satisfy the additional " free-fermion" conditions (FFC): 

f\u„0) = l, t = l,...,pN.2- (7) 

In this case, from @ and (|]) the corresponding eigenenergies of the Hamil- 
tonian (|ip for the case of free boundaries are given by 

PN-2 , . 

E = -2 ^ -cosr/ + cos— ^ , 1 < /tj < L - 1. (8) 
i=i ^ ^ ^ 

On the other hand, for periodic boundaries we have found solutions for the 
SU{3) model with r] = 2n/3, and for this case relations (|), (|) and give 
us 

Pi / 2nk \ 

E = - ^ 1 + 2 cos — , l<kj <L. (9) 

With the additional set of equations merged from the FFC (^, the num- 
ber of equations in and (0) exceeds the number of variables by pn-2- At 
the first glance we would expect no chance to satisfy the whole system given 
by @ and (0). But, surprisingly, the NBAE possess some hidden symmetry. 
In the next section, for the sake of simplicity, we consider the simplest case 
of the SU{3) model with q = exp(2z7r/3). 
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3. Free fermion spectrum for the SU{3) model with q = exp(2m/3): 

NBAE solutions. 

In 0, based on extensive numerical calculations we conjectured the ex- 
istence of free fermion solutions given by (§) and (^) for periodic and free 
boundary conditions, respectively. In the case of periodic boundaries these 
solutions happens in the sectors (no,ni,n2) = {k, k + 1, L — 2k — 1) (0 < 
/c < (L — l)/2) while in the open boundary case they belong to the sectors 
(fc, k, L — 2k) or {k, k + 1, L — 2k — l). In ^ exploring the functional equations 
related to the NBAE we were able to explain partially these conjectures in 
the periodic case. In the free boundary case we derived directly the NBAE. 
These solutions however do not belong to the numerically predicted sectors 
and do not satisfy the usual bounds po = no < Pi = Uq + rii < L. Moreover 
it is not clear if they correspond to non-zero norm eigenfunctions. In this 
section we present a direct explanation of these conjectures without the use 
of functional equations. 

The NBAEs for the anisotropic SU (3) Perk-Schultz model with anisotropy 
q = exp(2z7r/3), with roots satisfying the FFC, can be written as follows: 



PO 

n fi^^^^j) 


pi 

Y[ fivi,Uj) = 1, i = 1,2, . . . ,po, 


(10) 




i=i 




Pi 


Po 




n /(«i>«j\ 


)Ylf{ui,Vj) = l, z = l,2,...,pi. 


(11) 








f\u„0) = l, 


z = 1,2, 


(12) 



where po = uq, pi = uq + rii and we have used the relation 

F{x,y) = l/f{x,y), (13) 

which is valid for q = exp(^) or 77 = 
3a. Periodic boundary condition 

Before considering the general case let us restrict ourselves initially to the 
particular eigensector (1,2,L — 3). We have in this case pq = uq = 1 and 
Pi = Uq + rii = 3. For this simple case the first subsystem (^U|) consists of a 
single equation: 

f{v,Ui)f{v,U2)f{v,U3) = 1. (14) 
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From the definition of tlie function / given in (Q) we can show that the last 
equation (JT^) is equivalent to 



COs{v + U1-U2- Us) + COs{v - U1+U2 + U3) + 

+ cos(f — Ui — U2 + M3) = 0. (15) 

It is clear that this relation has the S3 symmetry under the permutation of 
the variables Ui,U2 and u^. Surprisingly it has also the 6*4 symmetry under 
the permutation of the variables Ui,U2,U3 and v\ 

The three equations of the second subsystem (|1T]) can be obtained from 
([T0| ) by the permutations v ^ Ui, so they also can be reduced to (pTSl) . Fixing 
the variables Ui, i = 1,2, 3, satisfying the FFC (O) and finding v from 



equation (^) we obtain the solution for the whole system (|lO]-0) in the 
eigensector (l,2,L-3). 

We show now that this method can be generalized to any sector {k, k + 
1, L — 2k — 1), supporting the conjecture 1 of our previous paper p. Let us 
consider the set {xj, z = 1, . . . ,po + Pi} of variables formed by the union of 
the two systems of variables: 

Xi = Vi (z = 1,2, . . . ,po), 

Xi+pg=Ui (i = 1,2, . . . (16) 

In terms of these variables the system of equations (|TD|-|I2p becomes 

r{xi) = f{xi,Xi), i = l,2,...,po+pi, (17) 
/^(a;i,0) = l, i=po + l,...,po+p,, (18) 

where r{x) = Ylf^^^ f{x,Xj). The first of these subsystems (|T^) possess 
Spo+pi permutation symmetry. We intend to show now that many of the 
equations in this subsystem are dependent, and for po = k and pi = 2k + 1 
we can also satisfy independently the FFC ([T8|). From (^) we may write 



_ sin(x - X, - n/3) _ b-qb, 
'^""^ }}, sin(x - x, + TT/S) - }}, qb - b, ' ^'^^ 

where for convenience we introduced the new variables 

b = exp(2ix), bj = exp{2ixj), q = exp(2i7r/3). (20) 
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The first subsystem (|T^ becomes now a system of algebraic equations: 

PO+Pl PO+'Pl 

Y[{bi-qb,)+ l[{qb,-b,)=0 z = l,2,...,po+Pi. (21) 

Using tlie standard symmetric functions: 
^0 = 1, 

Si = bi + b2 + ■ ■ ■ + bpg-^^p^, 



5*2 — 6l &2 + &1 ^3 + • • • ^po+Pi-l ^po+Pi' 



(22) 



Sm — ^ ^ii • • • ^im ; 

l<jl<i2<...«m<P0+Pl 



'S'po+pi — ^1 ^2 •••^po+pi! 

we can rewrite (^1]) as 

PO+Pl 

J2 (-l)"'^m&r°^'""'"(g" + gfo+Pi-™) = 0. (23) 

m=0 

Adding this last equation to the identity 

po+pi po+pi 
^-Po-Pi -Q _ ^^.) = ^ (-l)-S'„?^«+^'i-'"g-fo-m = 0, (24) 

j=l m=0 

we obtain 

PO+Pl 

J2 (-l)'"^m&?°+^'"™(g'" + + g-Po-Pi) = 0. (25) 

m=0 

For q = exp(2i7r/3) we have the following possibilities 

^ ^po+pi-m ^ ^-po-pi ^ 3g-po-Pi for po + Pi + m = 3n, 
g"" + gPo+Pi-™ + g-Po-Pi = for po + Pi + 7^ 3n. (26) 



Let po +Pi = 3fc + r, where k is an integer and r G {—1, 0, 1}. Inserting (^ 
into (|25| ) we obtain 

k 

^ (-1)^53 A"' = for r = 0, 

At=0 
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k—1 

J2hirS3,+2b;"' = for r = 1, (27) 

^J.=o 

k—l 

^(-l)^53^+i6r"^ = for r = -l, 

Ai=0 

where i = 1,2, ...,po + Pi- We then see that the subsystem (p!7|) can be 
satisfied if we impose 

S3u+p = 0, u = 0,1,..., k-1, (28) 

where p = for r = 0, p = 2 for r = 1 and p = 1 for r = —1. Since 5*0 = 1 
it is not possible to obtain solutions of this type for L = 3k {p = 0) and we 
have to limit ourselves to the cases where Po + Pi 7^ 3A;. 

Let us fix now pi variables Ui, satisfying the FFC (|TB|). Variables Vi, i = 
l,2,...,po can, in principle, be found from the system (p8D. In order to do 
that we use the decomposition of the symmetric functions depending on two 
set of variables, namely 

5*0 = 1, 

5*1 = Si + (Ti, 

82 = S2 + Si(Ti + (72, 

(29) 

min{m,pi} 
k=max{0,m—pQ} 
^PO+P\ ~ ^po'-^p\i 

where (Tj are the symmetric combination of the known variables Ui and Si 
are the symmetric combinations of the unknown variables Vj . Consequently 
the system (^) can be reduced to a linear system for symmetric functions 
Si, i = 1,2, ...,po. This system can be solved if the number of variables po is 
greater or equal to the number of equations k. Then we have the system 

Po = no, pi = no + ni, Uq < Ui, 

Po+Pi = 3A;±l, k<po. (30) 

These relations give us the constraint hq = k and rii = k+1, and consequently 
Po + Pi = + ni = 3k + 1, implying that solutions of (|28|) exist only for 
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p = 2. These solutions happen in the sector {k, k + 1, L — 2k ~ 1). Consider 
for illustration the case k = 2 {po = 2 and pi = 5). We have 5 variables 
ui,...,U5, which we fix with the FFC (|18|), and 2 unknown variables vi,V2- 
Using (p9|) we can write the system (^) as follows: 

82 = 82 + SiO"! +0-2 = 0, 

85 = S5 + S4(Ti + S3(J2 = 0. (31) 

The functions o'i,o'2,o'3,a4 and we know, so we have 2 linear equations 
for si = vi + V2 and S2 = fif2- 

3a. Free boundary condition 

The case of free boundary conditions is slightly more complicated. In- 
serting the definition into the subystem (p!7| ) we obtain the system 

po+pi 

sm{2xi - 27r/3) J] {cos{2xj) - cos(2xi - 27r/3)) + 
i=i 

Po+Pi 

+ sin(2xi + 27r/3) [] (cos(2a;j) - cos(2a;i + 27r/3)) = (32) 

i = 1,2, ...po+pi. 
As in the periodic case we may use the symmetric functions: 
^0 = 1, 

5*1 = cos2xi + cos 2x2 + cos2xp(,_|_pj, 

(33) 

S'po+p^ = cos2xi cos 2x2 ••• cos 2xpQ+pj, 

to rewrite ( P^ ) as 

PO+Pl / 

S'po+p,_„ sin(2x, - 27r/3) cos™(2xi - 27r/3)) + 

m=0 ^ 

+ sin(2xi + 27r/3) cos""(2xi + 27r/3))^ = 0, (34) 
z = 1,2, ...po+pi. 
On the other side we have the identities: 

PO+Pl 

sin2xj Y\_ (cos2xj — cos2xj) = 

PO+Pl 

= Spo+p^-rnsm2xicos^2xi, i = l,2,...po+pi. (35) 

m=0 
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Adding the equations ( P3| ) and (^) we obtain 



Spg+p^-m4>m{Xi) =0, i= 1,2, ...po + pi, (36) 

m=0 

where the functions 4>m{x) are given by 

1 

(pmix) = ^ sin(2xi + 27rp/3)cos"'(2xi + 27rp/3). (37) 

p=-i 

It is convenient now to find a reccursion relation for these functions. In order 
to do that let us consider a general sequence 

Vm = aa"^ + hp"^ + 07^". (38) 

This sequence satisfy the recurrence relation 

Vm+i + Ay^ + Bym-i + Cym-2 = 0, (39) 

if A, B and C are the coefficients of the third degree equation 

+ At^ + Bt + C ={t- a){t - /3)(t - 7), (40) 

with roots a, (3, and 7. Consequently our functions (t)m{x) satisfy the recur- 
rence relation: 

(t)m+l{x) + A(l)m{x) + 50„_i(x) + C(t)m-2{x) = 0, (41) 

where 

1 

t^ + At^ + Bt + C = n (^-cos(2x + 27rp/3)) = t^-3/4t-l/4cos6x. (42) 

p=-i 

Due to (|39|) one then have the recurrence relation for the functions (pm'- 

(i)m+i{x) = 3/4 0„_i(x) + 1/4 cos6x (f)m~2{x). (43) 

In order to iterate (^) we need also 4>o,(t>i and 02- Using the definitions (pTf ) 
we find 

00 (a;) = 0, 01 (a;) = 0, and 02 (x) = 3/4 sin6x. (44) 
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Let us give a list of several (piix), which we obtain by using the initial func- 
tions and the recurrence relations (^31): 



t>3{x) = 0, (p4:{x) = 9/16 sin6x, 4>rj{x) = 3/16 sin 6a; cos6x, 
t>G{x) = 27/64 sin6x, 4>7{x) = 9/32 sin 6a; cos 6a;, (45) 
hix) = sin6x (81/256 + 3/64 cos^6x). 



It is clear from ( PI) and (|I3| ) that 4>m{x) = sin6x$m(cos6x), where $m(^) is 
a polynomial. The degree fi of this polynomial depend upon the number m. 
One can prove by finite induction that for m = 3n + r, r = 1, 2, 3, 0m = if 
n — 2 + r < and otherwise it is a polynomial of degree /i = n — 2 + r. 

Returning to our main equation (^) and using the value we just obtained 
for n we easily verify that the left side of this equation is (up to common 
factor sin 6a;) a polynomial of cos6x with degree u = [(po+Pi ~2)/3], where 
[x] means the integer part of x. We can try to convert the equation (|36D 



into an identity on the variable x by equating to all coeffitients of the 
polynomial, i. e.: 

PO+Pl 

J2 Spo+pi-m $m(cos6x) = 0. (46) 

m=0 

We obtain k = u + 1 = [{po + pi + l)/3] linear equations on the symmetric 
combinations 5*^- Before considering the general case let us consider, for the 
beginning, the cases where po + Pi = 6, 7, and 8. Using ( ^5]) the relation 
among the functions (p and $ and the equation (H§), we obtain: 



5'4 + 3/4^2 + 1/4 5i cos6x + 9/16 = (po+Pi = 6), (47) 
^5 + 3/4 ^3 + 1/4 ^2 cos 6x + 9/16 Si + 

+3/8cos6x = (po + Pi = 7), (48) 
Se + 3/4 Si + 1/4 S3 cos 6x + 9/16 ^2 + 

+3/8 Si cos 6x + 27/64 + 1/16 cos2 6x = (po + Pi = 8). (49) 
We see that for Po + Pi = 6 and for Po + Pi = 7 one has two equations: 

^1 = 0, 

54 + 3/4 ^2 + 9/16 = 0, (50) 
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and 



S2 + 3/2 = 0, 

^5 + 3/4 53 + 9/16^1 = 0, 



(51) 



respectively. As far as po + Pi = 8 is concerned, the last term 1/16 cos^ 6a; 
makes impossible to convert ( 146|) into an identity. 

Let us consider the general case. In the same way as in the periodic 
boundary case we fix the variables Ui {i = l,...,pi) satisfying the FFC 
([T8| ) and try to find the variables f j (i = 1, . . . ,po) using k = [(po + Pi + 
l)/3] equations, which we obtain by imposing that all the coefficients of the 
polynomial (|46|) are zero. In order to find a solution the number of equations 
k can not exceed the number of variables po- The whole restrictions 



reduced to two variants: Uq = rii = k,pQ + pi = 3k and no = k, Ui = 
k + 1,pq + pi = 3A; + 1, which correspond to the sectors (fc, k,L — 2k) and 
(/c. A; + 1, L — 2/c — 1) respectively. One can show using the recurrence relation 
(|43| ) and the initial functions (|4^) that as long po +Pi 7^ 3fc + 2 we can obtain 
a consistent linear system for the symmetric functions Si . We use these linear 
equations in the same way as we have used (|2^) in the periodic case and this 
explains the conjecture 3 of p. For example, the systems ([50|) and (|5l|) give 
us free-fermion solutions for the sectors (2, 2, L — 4) and (2, 3, L — 5). 

4. Free fermion solutions of the NBAE for generic q 

In this section, differently from the results of last section, where the free 
fermion solutions were found for a specific value of g, we are going to present 
free fermion solutions of the NBAE that are valid for arbitrary values of q. 
These solutions are valid only for free boundary conditions and will happen 
in the special sectors of the SU{N)g model where the number of particles of 
each distinct specie is at most one, except for one of the species, that can 
be considered as the background (holes for example). The S{U)q symmetry 
ensures that the general sector containing all these solutions is the special 
sector where (no, . . . , n^v-i) = {I, . . . ,1, L — N + 1), that gives from the 
definition (|^) the values 



Po = no, pi = no + ni, no < ni, 
Po> [(po+Pi + l)/3], 



(52) 



Po = l,pi = 2, ...,Pn-2 = N - 1,Pn-i 



L. 



(53) 
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For this case the NBAE (|) consist of - 1 subsets (fc = 0,l,...,A^-l)of 
equations where the kth subset has precisely k + 1 equations. 

To illustrate our procedure let us consider for simplicity the SU (4) model 
in the sector (1, 1, 1, L — 3). In this case we have three subsets {k = 0, 1, 2) 
and three groups of variables {u^^\ u^^\ u^"^^): 

l = f{w,v,)fiw,V2), (54) 

F{vi,V2) = fivi,w) X fiVi,Ui)f{Vi,U2)f{Vi,Us), 

F{V2, Vi) = f{v2, W) X f{v2, Ui)f{v2, U2) f {V2, M3), (55) 

F(mi, U2)F{ui, Us) = f{Ui, Vi) f{Ui, V2), 

F{U2, Ui)F{u2, Us) = f{u2, f 1) f{u2, V2) , (56) 
F{U3, Ui)F{u3, U2) = f{u3, Vi) f{u3, V2) , 

where we write w instead of u'-^-', Vi instead of u^^^ {i = 1,2) and Ui instead of 

(2) 

u] (z = 1, 2, 3). Inserting (^ into the first equation ( p^ we obtain promptly 
two possibilities: 

sin2u7 = 0, (57) 
cos2t>i + cos2f2 — 2cosr7cos2it' = 0. (58) 

We do not intend to consider for the moment the first possibility (^Tf). 
Imposing the condition (|S^) we can check (eliminating, for example the vari- 
able w) that for {vi,V2,w) satisfying (|58D , besides (^4]) we have also two 
additional equalities, namely 

F{vi,V2) = fivi,w), 

F{v2,vi) = fiv2,w). (59) 
Consequently the second subsystem become: 

1 = f{vi,Ui) f{Vi,U2) fivuUs), 

1 = fiv2, Ml) /(f2, U2) fiv2, U3). (60) 

We will show bellow that the set of 5 equations formed by the second (|60D 
and the third subsets (0) contain only two independent equations, and con- 
sequently we may fix the variables Ui{i = 1,2,3) by imposing the FFC (|l^ 
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and find the two variables Vi and V2 from these equations. The remaining 
variable w is then obtained from (0). Consequently we find the free-fermion 
eigenspectra (|^) for arbitrary values of q or 77. 

The previous analysis can be extended for the sector (1,...,1,L — 
for general SU{N)q by exploring the general theorem: 

For any k fixed numbers Ui, U2, ■■■Uk one can find k—1 numbers Vi, V2, v^-i 
satisfying the two set of equations 

k 

f[f{vi,u,) = l {i = l,2,...,k-l), (61) 

i=i 

k k—1 

n F{u„Uj) = l[f{u,,Vj) (z = l,2,...,fc). (62) 

Leaving the proof of the above theorem for the moment, we can see that 
applying this theorem to {k — 1) known numbers fi, . . . ,ffc-i, one can find 
the numbers Wi, . . . , Wk-2 satisfying the equations 

k—1 

l[f{w,,v,) = l (2 = l,2,...,fc-2), (63) 

k—1 k—2 

n F{v„v,) = l[f{v„wj) {t = l,2,...,k-l). (64) 

Multiplying the ith {i = 1, . . . , k — 1) equation of the sets (^) and (0) we 
obtain literally one of the subset of the NBAE: 

k—1 k—2 k 

n F{v,,v,) = l[f{v,,w,)f[f{v„u,) {t = l,2,...,k). (65) 
Applying the above theorem k — 1 times we obtain a tower of numbers: 

Ul,U2,U3 Uk-2,Uk-l,Uk 

Vl,V2 Vk_2,Vk-l 

Wi, ...Wk-2 



yi,y2 

z 
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This imply that if we begin by fixing the k = N — 1 variables Ui {i = 1, . . . , k), 
satisfying the FFC (|T8|) we obtain by using recursively the construction (|65D 
the solution of the NBAE of the SU{N) Perk-Schultz model, with free bound- 
aries, in the sector (no, ni, . . . , un-i) = {1,1, . . . ,1, L — N + 1). The free- 
fermion like energies are given by (||) for arbitrary values of rj. The previous 
results (0) and (|60D are just consequences of the particular case where k = 3. 

Let us now prove the announced theorem. Let us fix {uj},j = 1, 2, k. 
The equation (|61|) can then be written as follows 

P{vi) = 0, {t = l,2,...,k) (66) 

where 

P{v) = Y\ ( cos 2uj — cos(2w — — n ( cos 2uj — cos(2t> + 7])]. (67) 



The use of the symmetric functions (^) (uj instead of Xi) allow us to write 
P{v) = j2 i-ir^'Sk-^ (cos'"(2^; - 77) - cos'"(2t; + r])) . (68) 

m=l ^ ^ 

Since cos(2t' — t]) — cos(2f + 77) = 2 sin 77 sin 2t; and a"^ — b"^ = {a — b) {a"^~^ + 
a™-26+ ... + 6™-i) we have 

P{v) = sin2t'p(cos2f ), (69) 

where p{t) is a polynomial of degree k — 1. We can factorize this polynomial, 
and apart form a multiplicative constant (f2)we can write 

fc-i 

P{v) = n sin 2v n (cos 2v - cos 2bi). (70) 

i=l 

Now ( |0U[ ) is easily solved: Vi = bi, i = 1,2, k — 1. 

Consider now the right side of the second equation (|62D . The relation 
( [70| ) allow us to write 



i=i j=i j=i cos(2ui + ?7) - cos26i 

P{2uj — rj) sin(2Mj + rj) 
sin(2uj — rf) P(2uj + rj) 



(71) 
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Using the expression ( |B7D for P{v) we obtain: 

^, U.i=i{cos2ui - cos{2uj - 2r])) - 11^^=1 (cos 2^^ - cos2nj) ^^^^ 

ni=i(cOs2Mj — COs2Mj) — nf=l(cOs2Mj — COs(2Uj + 27])) 

Since ni=i(cos2uj — cos2mj) = 0, we obtain 

n /(%. = -2i!f^i±4 n (73) 

sin(2uj — rj) fj-^ cos 2ui — cos(2uj + 27]) 



One can also check that 

sin(2Mj + 7]) cos 2uj — cos(2uj + 27]) 

sin(2Mj — 7]) cos 2uj — cos(2Mj — 27]) ' 

so that (|73D can be written as 



(74) 



k—1 k 

\{f{u,,v,)= n F{uj,u,), (75) 

1=1 i=l,i¥'j 

concluding the proof of the theorem. 

In Appendix B we derived the wave functions corresponding to the free 
fermion solutions obtained in this section. 

5. Free fermion spectrum for special values of q. NBAE solutions 
for the sector (1, 1, . . . , 1, 2, L - A^) 

For this sector the NBAE of the Hamiltonian (|1|) with free boundaries 
have about the same system of equations as in the previous section but the 
number of variables Ui is now k + 1 instead of k. We have in this case to use 
a modified theorem: 

For any k + 1 fixed numbers ui,U2, ■■■Uk, Uk+i one can find k — 1 numbers 
Vi,V2, ...,Vk-i SO that the two systems of equations 

l[f{v„u,) = l {i = l,2,...,k-l), (76) 

i=i 

k—1 

n F{u„u,) = llf{u,,Vj) {z = l,2,...,k), (77) 
17 



are satisfied for the special anisotropy rj = kTT/{k + 1). 

We begin the proof as in the case of previous theorem and as in (|69D the 
corresponding polynomial P{v) we obtain would have a degree k, however 
for the special anisotropy value rj = kn/ {k + 1) we can show that the degree 
of P{v) decreases to — 1, and we have the same situation as in the previous 
section. Consequently the proof follows straightforwardly in the same way as 
before and we have for the sector (1, . . . , 1, 2, L — N) the free fermion energies 
given by (§), with pn-2 = N and rj = {N — l)n/N. 

6. Free fermion spectrum for SU{N) model. Conjectures merged 

from numerical studies 

In this section we report some conjectures merged from extensive brute- 
force numerical diagonalizations of the quantum chain (|1]) with free boundary 
condition The analytical results presented in previous sections explain part 
of these numerical observations. In order to announce these conjectures let 
us define the special sectors 

S, = {[j;^],[^],...,[^^^],L~k), k = 0,l,...,L, (78) 

where, as before, [x] means the integer part of x. For example, for = 4 
and L = 7 the sectors are 

So = (0, 0, 0, 7), 5i = (0, 0, 1, 6), S2 = (0, 1, 1, 5), 
^3 = (1, 1, 1, 4), S, = (1, 1, 2, 3), ^5 = (1, 2, 2, 2), 
S'6 = (2,2,2,1), ^7 = (2,2,3,0), (79) 

and due to the quantum symmetry of the Hamiltonian we have the special 
ordering: 

5o C 5i C ^2 C ^3 C ^4 C 55 = ^6 D ^7. (80) 

This means that, for example, all eigenvalues found in the sector 5*3 can also 
be found in the sectors S'4, 5*5 and 5*6, and all eigenvalues appearing in the 
sector 5*7 can also be found in the sectors 5*6 and 5*5. The sectors 5*5 and 5*6 are 
totally equivalent. In this example let us call the sectors Sq, Si, S2, S3, S4, S^ 
as the left sectors and 5*6, S7 as the right ones. 

We can generalize this definition to any L = Nn + r, where n and r are 
natural numbers and r < N, obtaining L — n left sectors and n + 1 right 
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ones. Now we can formulate the conjectures merged form our bruteforce 
diagonalizations. 

CONJECTURE 1. For L = Nn + r (r = 1,2, N ~1) the Hamiltonian 
with p = and q = exp^in/N) has eigenvalues given by 

E, = -2 5:(cos(^)+cos(^)), (81) 

where I is an arbitrary subset of the set {1, 2, ...L — 1}. If k is the number of 
elements of the subset I and also Sk is a left sector then we find the eigen- 
values ^8l\ ) in the sectors Sk, Sk+i, ...Si-n- On the other hand if Sk is a right 



sector then we find the eigenvalues (pll) in the sectors S'l_„_i, Si^nj • • • , 'S'^+i. 



For L = Nn we have shghtly more dehcate picture. In this case we have 
the left sectors, the right sectors and a central one {n,n,n,n). 

CONJECTURE 1'. For L = Nn we can use conjecture 1 considering the 
central sector as a left or a right one. This is possible due to the coincidence, 
apart of degeneracies, of the eigenenergies in the eigensectors (n, n,n, . . . ,n) 
and {n — l,n,n, . . . ,n + 1). 

Consider now the a special subsets / = {1, 2, 3, . . . , k}, = 0, 1, . . . , L — 1. 
Due to the conjectures 1 and 1' the Hamiltonian (||) has the corresponding 
eigenvalues 

(h\ ■s^/ / 71" N /77I"nn , / TT , sin7r(2fc + 1)/2L , 
= cos +cos^ = l-2A;cos ^ — 82 

We can now formulate the remarkable conjecture: 

CONJECTURE 2. The eigenvalues ^8^ ) are the lowest eigenenergies in 
the special sectors ( [7^ j of the Hamiltonian with anisotropy q = exp{i{N — 
VjTi/N). Namely, for the left sectors we have Emin{Sk) = E^''^ while for the 
right sectors Emin{Sk) = E^''~^\ 
We can add here an additional conjecture: 

CONJECTURE 2'. The eigenvalues for k = L - n - 1 gives the 
ground state energy of the Hamiltonian: 

„ / ^ N / TT , sin7r(2n + 1)/2L 

E„ = l + 2(l-L + „)co=(-) ^—J^- (83) 

The last of the special sectors is 
Sl = ([^], [^], [ ^^^~^ 0), k = 0,l, L. (84) 
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In this sector we have only — 1 classes of particles and the Hamiltonian (|T[) 
is reduced effectively to a SUq{N — l)-invariant quantum spin model. The 
sector is of right type and due to conjecture 2 we can state the following 
conjecture: 

CONJECTURE 3. The eigenvalue 1 ^8^ ) for k = L — 1 gives the ground 
state energy of the SUq{N — 1) Hamiltonian with anisotropy q = exp{i{N — 
1)71 /N). This eigenvalue can be written as 

Eo = -2(L-l)cos(^). (85) 

For = 3, for example, we get the XXZ model with anisotropy A = —1/2 
and Eq = 1 — L, a result that was first observed in and produced quite 
interesting consequences 0, ^ . 

7. Summary and conclusions 

Most of the researches related to exact integrable systems ends with the 
derivation of the Bethe-ansatz equations, whose solutions provide, in princi- 
ple, the eigenvalues and eigenvectors of the Hamiltonian or transfer matrix 
of the associated model. These equations are highly non linear and in several 
cases thanks to some appropriate guessing on the topology of the roots in the 
complex plane, they are transformed in the bulk limit (L —>■ oo) into inte- 
gral equations, that provide most of the thermodynamic physical quantities, 
whose calculations depend upon eigenvalues only. Even in this limit theses 
simplifications does not help for the calculation of quantities depending di- 
rectly on the eigenvectors like form factors and correlation functions, since 
the amplitudes of the eigenvectors, obtained in the Bethe basis, are given as 
combinatorial sums of plane waves that are highly nontrivial to operate in 
the thermodynamic limit. 

Since the exact integrabihty is a lattice size independent property, the 
exact solution of the associated Bethe-ansatz equations is a quite important 
step towards the complete mathematical and physical understanding of these 
models. Due to the complexity of these equations, the solutions are known 
analytically in very few particular cases, like the free fermion case and the 
XXZ chain with anisotropy A = —1/2 0. It is interesting to mention that 
in this last case only some ratio among the amplitudes of the eigenvectors 
are known and these results provide a quite interesting connection of this 
problem with that of the sign-alternating matrices 0]. In a previous paper 
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we presented some new solutions of the associated Bethe ansatz equations 
of the simplest generalization of the XXZ chain with two conserved global 
quantities, namely the SU{3)q Perk-Schultz model with periodic and free 
boundary conditions. In this case the Bethe ansatz equations are of nested 
type (NBAE) and these solutions were obtained for the special anisotropy 
value q = exp(z27r/3), being the first example of analytical solutions of NBAE 
for finite lattices. The results in were derived by exploring the functional 
relations merged from the associated NBAE, and presented two difficulties: 
the solutions for the open lattice happen in sectors out of the usual bounds 
(no < ni < n2), being not clear if they do not correspond to zero- norm 
states, and the precise sector where the periodic solutions happen is not 
known. In §3 of the present paper, in a different way, by adding to the set of 
NBAE the FFC (|I^) we were able to solve these difficulties by showing the 
existence of free-fermion solutions in the sectors satisfying the usual bounds 
(no < ni < 77-2) for the periodic and open cases. On the contrary to the 
functional method used in [§, the success of the direct method presented 
in §3 for the SU{3)g case enable us to generalize these solutions for some 
sectors of the Perk-Schultz SU {N)q models, revealing the existence of free- 
fermion branches, at the special anisotropy q = exp(z27r(iV — 1)/A^). It is 
important to mention that our numerical calculations reveal other energies 
beyond those predicted in the calculated free-fermion branches. Although 
we did not prove analytically, an extensive numerical analysis shows that 
the low lying energies, including the ground-state energy belong to these 
classes of free-fermion like solutions we have found. In fact it was these 
numerical results that induce our analytical calculations. The numerical 
analysis indicated a set of conjectures presented in §6, part of them were 
proved in §3 - §5 and appendix A. 

After the success of finding free-fermion solutions at special values of 
anisotropics of the SU{N)q model, we tried to find other solutions for other 
anisotropics, at least for some simple sectors. Surprisingly we were able to 
find free-fermion solutions for all values of the anisotropy in the sector were 
we have at most one particle of each specie, except for one of the species (the 
background one) . In §4 we show the existence of these solutions for the most 
general sector within this class, i. e., (no, . . . , n7v_i) = (1, . . . , 1, L — -|- 1). 
Beyond the free-fermion energies degenerated with the lower sectors, appear 
in this sector the energies given by the sum of pn-2 = N — 1 free-fermion 
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energies ej, i. e., 



N-l 



N-1 



Xj 



where 



Xj — exp(i— 1 < kj < L 
Ij 



1. 



(86) 



(87) 



These eigenvalues expressions, vahd for arbitrary values of g, induce us to ex- 
pect that a direct analytical derivation of the related eigenfunctions, without 
the use of the NBAE might be possible. In fact in the appendix B we derive 
these eigenvectors, which are given by combinations of Slater fermionic deter- 
minants of one particle wave functions \E'j(m;). Denoting the configurations 



where the Zth particle is at position mi (l — 1, . . . , N 
I mi, . . . jTUn-i > the eigenfunctions are given by 



l;mi — 1, . . . , L) as 



\'^{X1,...,XN-1} >= J2 (1 

mi,...,mjv-i=l 

^i(mi) ^i{m2) 



-/(mi,...,mjv_i) 



xdet 



^i(miv_i) 
*Ar-i(mAr-i) 



|mi, . . . ,mAr-i >, 



where 



^ .(m) = 1 - — - (1 - qxj) jx 



X 



(88) 



The factor /(mi, . . . ,m7v_i) is an integer given by the minimum number of 
pair permutations of the set (mi, . . . ,mAr__i) (m'^, . . . ,m^_^), such that 
rn!^ < m'2 < ■ ■ ■ < m'j^_i. It is interesting to observe that at the isotropic 
point the wave function amplitudes are just given by standard free fermion 
Slater determinants. 

In §5, continuing our search for free-fermion like solutions we were able 
again to find solutions in the sector (no, . . . , nM-i) = (1, . . . , 1, 2, L — N). In 
this case the free fermion solutions happens only at the special value of the 
anisotropy r] = {N — l)n/N. Although we believe that should be possible to 
obtain free-fermion like solutions for arbitrary eigensectors, this procedure is 
not straightforward, and each case deserving some ingeniousness. As another 
example we show in appendix A the existence of free fermion solutions in the 



22 



sector (no, rii, n2, n^) = (1, 2,2, L — 5) of the S'f/(4)q model at the anisotropy 
q = exp(237r/4). 

Ahhough we have only derived wave functions for the sector (no, . . . , nAr_i) 
= {1, . . . , L — N + 1) we expect to obtain in the future the corresponding 
eigenfunctions for the other sectors. This program is very important, since for 
example, it may teach us to derive in a closed form the wave function for the 
interesting problem of the A = — 1/2 XXZ chain, or the sector (0, n,L — n) 
(n = 0, . . . , L) of the SU (3)g model at g = exp(i27r/3). 
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Appendix A. SU(4) model with q = exp(3i7r/4), 
sector (1,2,2,L - 5). 

The nested Bethe ansatz for SU (4) Perk-Schultz model with free bound- 
ary can be written as follows (see (0) and (|^)): 



PO Pl 

Yl F{wi,Wj) = Y[f{wi,Vj), i = 1,2, ...,po, (89) 

Pl PO P2 

P2 Pl 

n F{u„u,) = l[f{u„v,), z = l,2,...,p2 (91) 
/^(n„0) = l, t = l,2,...,p2, (92) 

where we restricted ourselves to the solutions satisfying the FFC (|^). 

We have found free-fermion solutions for arbitrary values of q in the sector 
(no, ni, n2, n^) = (1, 1, 1, L — 3) (see §4) and for q = exp(3z7r/4) in the sector 
(1, 1,2, L — 4) (see §5). We intend now to obtain solutions in the next simple 
sector (1, 2, 2, L — 5). For this case po = 1, pi = 3, p2 = 5. Choosing the five 
variables Ui satisfying the FFC (^2]), we look for three variables Vi and one 
variable w, satisfying the nine equations given in (|89| ) - (|9ll) . 
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Let us begin with the third group of equations (^Tj) . The left side of these 
equations can be written as L{ui), i = 1, . . . , 5, where 

iM=^n/(«.«i) (93) 

Using the definition (|^) and fixing the anisotropy parameter rj = Sn/A we 
obtain after tedious but straightforward calculations: 

L(u) — If \ 2/ 

S S5 + S3 + Si - C(S4 + S2 + 1) - C (S3 + 2si) + 



L{u) + 1 

+C^(S2 + 2) + C^i - I S4 + S2 + 1 - C(S5 + S3 + Si) - (94) 



-1 

-C^(S4 + 2S2 + 3) + C%S3 + 2si) + C^(S2 + 3) - C^Si - 

where we have used short notations for the symmetric functions 

si = cos 2ui + cos 2^2 + ■ • ■ + cos 2-U5 

S2 = cos 2ui cos 2^2 + ■ ■ ■ + cos 2^4 cos 2m5 

(95) 

S5 = cos 2ui cos 2m2 ■ ■ ■ cos 2^5 

and for the trigonometric functions: c = cos2u and s = sin 2m. It is clear 
that 

5 

(cos 2uj — COS 2u) = S5 — c S4 + S3 — S2 + Si — = 0, (96) 



for u = Ui. Using the last identity we can exclude S5 from (plf). After some 
simplifications we obtain the more simple fraction: 

L{u) - 1 _ S3 + Si - C(S2 + 1) _ . - 

L{U) + 1 %4 + ,2 + l-c(s3 + Si)-c2 + c4' «l---«5.iyO 

On the other hand the right side of the third group of equations ( |9TD can 
be written as R{ui), i = 1, . . . ,5, where 

Riu) = n/(«,^.)- (98) 
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Again, rather straightforwardly we obtain: 



R{u) - 1 ^ 272^2 + 4c ai + 72(1 + 2c^) 

R{u) + 1 ^4a3 + 2(Ti + 2v/2c(T2 + \/2(3c-2c3)' 

where we have used the short notations for the symmetric functions 
CTi = cos2t>i + cos2f2 + cos2f3, 

cr2 = cos 2t>i cos 2t>2 + cos 2t>i cos 2vs + cos 2^2 cos 2vs, (100) 

(73 = COS 2t>i COS 2t>2 COS 2t>3. 



By equating the right sides of equations (pTl) and (^) we obtain a polyno- 
mial equation of degree 6 with respect to the variable c. Using ( p6|) we can 
exclude the terms with and and obtain a polynomial equation of the 
4th degree with respect to the variable c. It is amusing to verify that all the 
five coefficients of this polynomial become zero if we choose 

_ 2Si + S3 + Si S2 

+ (101) 

/ i + S2 

S3 + 2s5 — Si S2 — 2si S4 
~ 2v^(l + S2) ■ 

This result means that for arbitrary values of Mi, . . . , M5 we calculate by using 
(p5| ) si,...,S5, and the above relations ( p.01| ) together with the definition 
(|100|) enable us to find the variables f 1, f 2, as the roots of the equation 

COS^ 2t> — CTi COS^ 2f + (72 cos 2^; — (73 = {v = Vi,V2,V3). (102) 

Let us consider now the first equation: 

3 

Y[f{w,Vj)-l = R{w)-l = 0, (103) 

where R is given by (^). Due to ( pOj) the equation ( |103| ) becomes: 

2^2(72 + 4 cos 2w cTi + V2{1 + 2 cos^ 2w) = 0. (104) 
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Inserting ( [I01| ) into this last equation we obtain 

(si — cos 2w) ( ^iJl^ _ cos 2w\ = 0. 
^ V 1 + S2 / 

We have two possible solutions for cos 2w. Choosing the solution 

Sl + S3 



cos2w 



I + S2 



(105) 



(106) 



we can check by using standard algebra that for this choice of w and Vi,V2, 
obtained from ( |102|) all equations of second group (^) are satisfied automat- 
ically! 

As a result, we explained one more free-fermion branch found in our 
numerical observations. 



Appendix B. Eigenvectors for the eigenvector in the sector 

(1,...,1,L-Ar+1) 

In §4 we found a set of solutions in the sector (1, . . . , 1, L — + 1) of 
the SU{N)q model with an arbitrary value of the anisotropy q, and free 
boundary conditions. The corresponding eigenenergies are given by (||), and 
can be written as: 



N~l N~l 1 

j=i j=i y 



1^ 
J 

Xj 



(107) 



where 



Xj = exp[i 



irkj , 



l<kj<L-l. 



(108) 



We intend to prove that the corresponding eigenvectors have the following 
form 



\i^{xu-,XN-i} >- 



L 

mi,...,miv-i=l 



-/(mi,...,mjv_i) 



xdet 



^i(mi) ^i(m2) 

^2("^l) ^2("^2) 
^'Ar_l(mi) *Ar_i(m2) 



N-1 



{rriN- 



mi, . . .,mN-i >, 
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where |mi, . . . ,m7v-i > are the vector basis representing the configuration 
where the ith particle is located at site rrii {i = 1, . . . , N — 1, 1 < rrii < L) . 



The Slater determinants entering in ( p.09| ) depend upon the set of one-particle 
amplitudes \E'j(m), that apart from a normalization factor are given by 

^j{m) = - xf - (1 - qxj) /xf. (109) 

The factor /(mi, . . . ,m7v-i) is an integer number that depends on the rel- 
ative order of the sequence {mi,m2, . . . ,mAr_i}. We are going to find this 
dependence in the procedure of proving ( |109| ). 

The application of the Hamiltonian (|I]) to the above vector leads to a set 
of equations for the amplitudes \Ef(mi, . . . , mAr_i). 

Consider initially the amplitudes where \mi — mj\ > 2 for all pairs 
These amplitudes give us a set of equations ("regular" ones) for one-particle 
amplitudes that are satisfied for every product in the determinant separately. 
Consider now the case where rrii = m and rrij = m + 1 for precisely a single 
pair {i < j). The equation in this case are given by 

N-l 

E^{...,m,...,m + 1,...) = - [^(•••,"^fc-l,---) 

+^(. . . , mfc + 1, . . .)] - ^(. . . , mi, . . . , m + 1, . . .) 
— \E'(. ..,m...,m + 2) — \E'(...,m + l,...,m, ...) 
+ [{N - l)q +{N- 2)/q] ^(. . . , m, . . . , m + 1, . . .)• (HO) 

Comparing this equation with the "regular" ones, derived previously we ob- 
tain 

\E'(. . . , m, . . . , m, . . .) + \E'(. . . , m + 1, . . . , m 1, . . .) 

-^(. . . ,m+ 1, . . . ,m, . . .) - i^(. . . ,m, . . . ,m + 1, . . .) = (111) 



Inserting (|109|) into ( |1 1 1| ) we verify that the first two terms vanishes while 
the third and fourth ones correspond to the same determinant, apart from a 
minus sign. As a result we obtain an equation for the factor / in ( |109| ): 



/(..., m, m + 1, ...) = -1 + /(..., m + 1, m, .. .). (112) 



27 



We are now ready to fix the function /. Let f {mi, 1712, . . . ,m7v-i) 
all {rrij} satisfying to the order: 

nil < 1712 <■■■ < ^N-l- 

Then ( |112| ) fixes / for all remaining configurations. 

Consider, for example, the SU{4:)q case. We have 6 possibilities 



0, for 



;ii3) 



/(mi < 1712 < ITT-s) 

f{m2 < rrii < m^) 
f{m2 < m-i < mi) 
/(ma < m2 < mi) 





/(mi < m3 < m2) 
/(ms < mi < m2) 
3. 



It is clear that / is equal to the minimal number of pair transpositions nec- 
essary to put the configuration in the order mi < m2 < ms. 

One can easily check that the remaining equations coming from the other 
amplitudes reduce to the already obtained relation ( |112| ), concluding the 
proof that (|109| ) are the eigenfunctions corresponding to the NBAE solutions 
derived in §4, with eigenvalues given by (|107|). 
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